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ABSTRACT 

Solar stills similar to the type first built in Ias Salinas, 
Chile, in 1872, are today the type which show the greatest promise 
for utilization of solar energy for the conversion of sea water to 
fresh water, 

In this paper, a method of improving the efficiency of such a 
solar still by cooling the condensing surfaces with sea water is 
investigated. Results are negative in that by so doing the efficiency 
of the still is decreased, 

A dimensional analysis provided a possible correlation of the yleld 
of a solar still with the temperature of the water in the evaporating 
tray, and see temperature of the external glass surface. This correla=- 
tion should prove useful in future efforts to design and operate solar 
stills at maximum efficiency, and might have applications in related 
fields of heat transfer, 

The writer, with deep appreciation, wishes’to thank Professors 
P. F, Pucci and C,. P, Howard for their invaluable advice and counsel, 
and Mr, R, P, Kennicott for his assistance in the construction of the 


solar stills used in the investigation. 
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LIST OF SYMBOLS 


A Area, rt”, also indicates Still "A" 
B oti.) eB 

Dy Diffusivity £¢°/Hr 

g Acceleration, 4.17 »° 10° £t/Hr* 

h Coeff, of heat transfer pru/arrt“* . 
K Thermal conductivity BTU/Hrft oF. 
L Characteristic length; 3.915 feet 

fi Mass flow rate, lbs/Hr. 


Di fonl ~ GO . -R) . PRCpD 
N mensionless parameter ale Ts) (Fr-R) —~— 


Ty ae 
Nt Lumped constants and temperature variables of N, = Gp(12)? ed 
A Ple Lewis 
P Vapor pressure (see subscripts) 


Pe Peclet number, nC, = 4,15 m 
Ki 


Q Quantity of heat, BIU 
Rate of heat flow BTU/Hr 
7 Temperature, (see subscripts) 
— some characteristic angle, radians 
B Coefficient of expansion 1/T 
€ Emissivity 
A Latent heat of evaporation, BIU/1b. 


Viscosity lbs/Hrft. 


mS 


Density, 1bs/ft? 
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LIST OF SYMBOLS (continued) 
Subscripts 
a Ambient and average. Ta indicates ambient temperature 


r Film. Tp indicates film temperature = Ty + T, 
2 


g Glass, external, fT... indicates an area-weighted 
average external glass temperature for a still, 
P, indicates saturated vapor pressure at De 


n North : Tn indicates temperature of external. 

Ss South : glass surface of North side of still. 

e Kast : Pn indicates the saturated vapor pressure 

W West : at Tn, A, indicates the North face of 
Still A. 

t Tray. Ty, indicates temperature of water in the 


evaporating tray. P,; indicates saturated vapor 
pressure at qT. 
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1. Introduction, 

Almost from the time when it was considered a god, the sun has 
challenged man to use its boundless energy to perform some of man's 
many tasks, Scientists and homeinventors alike have produced endless 
devices to answer this challenge, but rarely have they become more 
than curiosities, 

A similar challenge has come from the sea, That a desert, a 
potential garden land but for lack of fresh water, can exist by the 
side of a vast body of sea water is still a problem of world=wide 
ampOrt. 

In 1872 man answered one of these challenges by means of the 
other. In Ias Salinas, Chile, a solar sea~water distillation plant 
of over 6000 gallons a day yield wae constructed, and it operated 
successfully for many years. Since that year countless varieties of 
solar stills have been suggested or built. However the Chilean-type 
Btills, with certain improvements, maintain their position today as 
the type of solar still offering the best yleld-to-cost ratio. (25) 

The economics of solar energy present a paradox in that solar 
energy is "free", and the means are at hand to collect it, yet a 
solar solution to an energy problem is often the most expensive solu~ 
tion, Such is the case in the United States. Therefore, although 
there are indications that the sea will (and possibly must) begin to 
supply some of our fresh-water needs, the energy source will probably 
not be solar, The reason for this is twofold, First, the cost of 
the equipment per gallon a day collected, and secondly, the amount of 
land area required. It is estimated by lof (25) that the lowest 
ultimate cost per 1000 gallons of solar water in this country will not 
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be less than one dollar, which compares poorly with the present five 
cents to twenty cents per 1000 gallons from municipal sources, The 
average intensity of solar radiation in this country is 1500 BTU/day, 
£t* (25). If a solar still operates at 60 percent efficiency, about 
ten square feet of collector area are required for each gallon per day 
produced, A community area the size of the Monterey Peninsula would 
then require roughly 1800 acres of stills to supply its needs, If 
multiple effect distilzation were employed, this last figure could be 
quartered, but the cost per 1000 gallons would be trebled. (26) 
Despite the above conjectures on the future of solar distillation 
in our own country, in other countries where labor conditions, material 
and land availability are far different than in our own, a solar dis- 


tillation is looked upon as a method holding great promise (9). 


e les Salinas SAS 


The Ias Salinas<type still can be most easily described as 
resembling a long pup tent, with glass ends and roof, oriented so 
that the ridgepole is pointed Hast and West. A thin layor of sea- 
water is held in a dull black tray - the floor of the tent. 

Glass has unique properties that specially recommend it to uses 
involving solar heat collection, I1t is almost completely transparent 
to radiation of the solar wave lengths, absorbing, dependent on impurity 
content, as low as one percent of the total incident energy, Reflection 
losses are higher, being about eight percent of the incident energy for 
@® beam normal to the glass, However, glass is virtually opaque to radia- 
tion from surfaces below 300° to 400° F,, and in this region has an 
emissivity of about .96. 

Therefore, depending on incidence angle, up to 91 percent of the 
sun's energy is transmitted to the collector tray of the Ias Salinas 
still, which can reach a temperature of 140° ~ 180° at midday, Re- 
radiation from the tray is determined only by emissivities and the 
temperatures of the water in the tray and the glass roofs, This is 
about ten percent of the incident energy. Water, evaporating from the 
tray, condenses on the sloping glass roofs (which are generally 25° - 
30° F, cooler than the tray at midday), and is collected, In cooling 
these condensing surfaces, radiation to the sky and natural convection 
are equally important; wind is more important than either, This perhaps 
can be best illustrated by the heat balance included in appendix (1). 
Although the amount of incident radiation was not measured, it is 
believed the figures are approximately correct, 


Maximum yield that has been reported for such a still is .18 gallons 
5) 





a day per square foot of tray surface (Algeria), Tuxis varies with 
latitude and ambient conditions such that the maximum reported yield 
in this vicinity is one-tenth gallon per day per square foot of tray 


surface (Berkeley), 





2e Squipment and Instrumentation, 
Construction, 

The solar stills used in this investigation are of the Ias Salinas 
type. Since no method of measuring the incident energy was at hand, 
two stills were constructed of exactly equal description, one to be used 
as a standard (henceforth called Still "B"), measuring the effect of 
alterations to the other (henceforth called Still "A"),. They have a 
tray dimension of 47" x 47", The glass roofs and ends are three=six- 
teenths-inch window glass, the roofs sloping at an angle of 45°, The 
bottom of the still was covered with a sheet of one-eighth=inch press- 
board, A twor~and-onewmhalf-inch thick fiber glass blanket lies between 
this bottom and the tray. The frame of the stills is redwood, and all 
joints, both wood=wood and wood-glass were sealed with Pabco Hydroseal, 
The galvanized sheetmetal tray was painted with Sherwin Williams dull 
black Enameloid, 

Figure 1 shows the general appearance, and indicates the more 
important dimensions of these stills. Figure 2 is a photograph of the 
actual installation, 

The water level in each still was adjusted so that the evaporating 
tray held 40 lbs. of water. This was the minimum amount that would 
prevent some portion of the trays from being exposed toward the end of 
the day, The depth of water varied, due to tray buckling, from one- 
fourth to three-fourths inches, Condensate was collected from each 
condensing surface separately. Collecting troughs, milled into the 
redwood frame of the stills below the condensing surfaces, carried the 
condensate to external collecting bottles, 

Since data was tote comparative, no attempt was made to design 
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the stills for maximum efficiency, Available glass with a greenish 
tinge (indicating presence of iron) was used, Likewise it is possible 
that & thinner glass could have been used with some efficiency increase, 
The length to width ratio was not as large as desirable, but this per- 
haps increased the overall efficiency due to a better angle of incidence 
in the afternoon, 

Still "A" was modified to permit cooling the North and South con~ 
densing surfaces with water, 

Two tanks of combined capacity of 70 gallons were used as insulat= 
ed storage tanks. One sheetmetal tank of equal capacity was used as 
a receiving tank for the cooling ater discharge and for the water 
used to flush the stills, Painted black, it was also used as a radia- 
tor to cool the cooling water each night. In the morning the water, 
ata temperature of 48° - 55° F., was drained into the insulated tanks 
where it was available for further cooling purposes, Pressure on the 
cooling system was maintained by a small air compressor, 

Water from the pressurized tanks was regulated by needle valve and 
gent to Still "A® by insulated one-fourth inch Saran Tubing. At the 
still it was sprayed over one or both of the sloping faces from three- 
eighths-inch Saran Tubes fixed across the top of the North and South 
faces, Several methods of evenly distributing the water were tried, 
Small holes spaced evenly the length of the tubing were unsatisfactory, 
due to clogging. The method finally used was to slit the plastic on 
one side the full length of the glass and then insert a felt gtrip two 
inches wide into the slit, With this method the system was non~cloge- 
ing and could be evenly spread across the top of the glass, However, 


the flow tended to concentrate into streams, and at about midway down 
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the glass only about 25 percent of the glass was covered by the water, 
Had the method been successful in increasing yield, a better means 
would have to be found to distribute the cooling water, 

There were also filling and discharge lines leading to each still 
to enable flushing, Both flushing and cooling could be regulated to 
occur at any desired tray temperature by means of Solenoid Valves 


operated by a thermostat beneath the tray of Still "A®*", 


Comments on Construction, 

(1) It is felt that Hydroseal is not an ideal sealant since it imparts 
& strong, unpleasant taste and aroma to the water. 

(2) It is felt that a still can easily be designed which will minimize 
the need for sealant by recessing the glass edges in grooves, and by 
careful attention to the design of the lower edge joint. Figure (3) 
shows the type of joint used in subject stills, which is considered to 
be improper; and shows a better type of jointure which would minimize 


leakage and the contact of the condensate with the sealant, 


location of Stills, 

The stills were located on the roof of a three-story building at 
the U. S, Naval Postgraduate School, Monterey, California. Trees in 
the area to the North and East provided a wind break so that during 
the investigations described the local wind was rarely above seven 
knots, The stills, however, had a clear view of the sky in all direc- 
tiong with one important and awkward exception: one tree was so located 
that after 15 March Still "B" did not clearly see the sun until 0820, 


and Still “A® untill 0840, 
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The trays were about 15" above the rooftop. Still "A" was ten 
feet Hast of Still "B", Both were oriented so that the sloping roofs 


were facing North and South, 


Instrumentation, 

Both stills were equipped with copper~constantin thermocouples to 
enable recording of the following temperatures: 

Tray temperature; thermocouple soldered to center of tray. 

Water temperature; three-sixteenths-~inch from water surface, 
near tray center and where water was one> 
half inch deep, 

Vapor temperature; inside aluminum radiation shield at west 
end of stills, midway between ridge and 
traye 

Glass temperature, exterior; centered on the North and 
south faces and centered on the South 
triangle of each Hast and West face, 

Cooling water ("In" and "Out"); Still "a", 

Ambient air temperature, 

Even though thermocouple terminal boards were at similar locations 
on the West end of each still, and shielded from the sun, it was found 
that temperature of these boards varied as much as four degrees; and 
that initial use of only one ice junction, at Still "B", caused tempera- 
ture errors of the same magnitude due to secondary thermocouple effect 
at these boards, 

Three terminal boards were required at each still and eventually 


each therminal board was connected to an individual ice junction, 
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Temperatures were recorded upon a Honeywell Multiple Point Recorder 
using a constant~balance type of circuit. 

Millivolt readings could be accurately taken to the second decimal 
place, with the third place estimated, The recorder was located ina 
third=-floor room and connected to the thermocouple terminal boards by 
30-foot copper leads, 

Wet and dry bulb temperatures were taken hourly during test runs 
with a sling psychrometer, Wind readings were similarly taken ata 
point between the two stills using a hand-held anemometer, U, 5. Navy, 


type AN/PMQ=-3, (Made by Bendix Corp.) 


Comments upon Instrumentation: 
(1) Additional foresight would have dictated placing thermocouple ter= 
minal boards on North side of stills where shielding from the sun would 
have been less difficult. 
(2) Unfortunately the temperature information cannot be viewed with 
complete confidence, since toward completion of the investigations 
described, evidence was found of ground looping between the tray and 
water thermocouples, Maximum errors occured at midday, and it is 
believed that in some cases they approached two degrees, Final runs 
were made with the tray thermocouple disconnected, giving better re~ 
sults, 
(3) Since several investigators have reported difficulty in finding 
a method of attaching thermocouples to the glass surfaces of the still, 
the following description of the method used is included: 

After numerous failures involving Du Pont and Glyptal cements, 


Glyptal varnish, Armstrong epoxy resin cements and several plastics a 
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method was adopted that produced a lasting and what is believed to 

be firm connection between thermocouple bead and the external glass 
surface, A small cube of rubber weatherstripping, one~fourth inch on 

a side, was coated with Goodyear Pliobond to waterproof it. Then a 
small piece of aluminum foil about five-thirty~seconds inch square 

was centered on one face of the cube, Pliobond was applied to the 
periphery of this face and to a spot on the glass, leaving a bare cen- 
ter on each where the thermocouple bead was to sit. With the bead in 
the glue-free center, the cube was pressed onto the glass and taped 
down to provide pressure while the glue dried, Once dry, the tape 

was removed and excess glue was scraped off the glass. This type of 
installation has withstood many rains and nightly dews without loosen=- 
ing, Not only is the bead held firmly against the glass by the elasti- 
city of the rubber, but the aluminum square firmly hues both the bead 
and the glass, giving an effectively increased contact area, The ther- 
mocouple is also shielded from the sun's radiation, It was initially 
assumed that, although the rubber cube insulates the contact area of 
the glass from the air, temperature equalization occurs to the sides 

to a satisfactory degree. A subsequent examination by relaxation plot 
ting (with conservative assumptions) showed that for a solar input of 
265 Bru/ft* hr, the thermocouple should read about four-tenths of a 
degree higher than the actual glass temperature, and the discrepancy 

ig believed to be actually much less than this, 

(4) The assumption that one thermocouple placed at the center of a 
glass face would give an acceptable mean temperature for that face is 


considered to be valid, since measurements with a thermocouple similar 
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to that described in 3 above, but mounted on a hand probe, showed 
variation across a glass face in the horizontal and vertical direction 


to be in the order of two-tenths of a degree, 
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3. Experimental Procedure, General, 

The investigation can be discussed under three interrelated but 

separate topics: 

a) Observations to determine relative performance of each 
condensing surface in respect to yleld, and the perform- 
ance of Still "A" relative to Still "B" under normal operat- 
ing conditions, 

b) Observations to determine the effect on the yield of Still 
"A" of changing the heat transfer conditions of the con- 
densing surfaces, 

c) Observations to determine the relation of the yield of 
the condensing surfaces to the dependent variables of 
tray water temperature, and external temperature of the 
glass condensing surfaces, 

Certain experimental procedures were common to each of the above 

phases: 

a) Until 15 March, the stills were flushed automatically during 
the night, and the measurements of yield were made on the 
basis of a 24-hour run, However, by March 15, the sun had 
moved northward enough to cause the aforementioned tree 
shadow to fall across the stills, Still "A" did not clear- 
ly see the sun until 0820, and Still "B® until 0840, This 
unavoidable shielding necessitated starting comparison 
tests after this date at 0900, and made exact duplication 
of initial conditions in each still difficult. This was 
achieved as closely as possible by a combination of flush=- 


ing both stills and covering Still "B" with a tarpaulin 


are 





v) 


c) 


between 0820 and 0840, By May 13 the shadow crossed the 
face of the stills at approximately the same time and for 
an equal period, at 0730 Pst and no equalizing procedure 
was required on that date. 

Ail daily yield measurements were made with a spring scale 
at 0800. 

On days when the mass flow rate was determined, yield 
measurements were made at 30-minute intervals, using a 
500 or a 1OOML graduate, These measurements were con= 
verted to mass flow rate, using central differencing 
methods, by means of the following numerical differentia- 


tion forma (3), 


L2hD = Yfi-z) — 8YE-)) + B8Yhs) — You) 


For easy application this formula was converted to 


~Vu3g9 + 7% + WKiryr Vlireo 


12hm 


where? i = Mass flow rate, lbs hr at time 


h = interval between measurements = 1/2 


Vie) 


the lst backward difference, = amount (lbs) 


u 


of condensate produced during the 30-minute 
period preceding time(A) , similarlyYiuieg = 
amount of condensate produced during the 
30—minute period preceding time (< + 30). 


This formula should produce errors of the order n+, 


ad) Temperatures were also recorded at 30-minute intervals. 


e) For convenience, fresh water was used both for cooling 
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water and for the distilling trays. 


It is considered that the results would not be affected by this 














4, Observation to determine relative performance of each condensing 
surface, in respect to yield, and the performance of Still "A" 
relative to Still "B" under normal conditions of operation, 

During the period of the investigation, yield measurements were 
taken with both stills operating under as nearly identical comitions 
as possible, These measurements were used as a basis for comparison 
when alterations were made to the characteristics of Still "A", and 
also to determine the relative capability of the various condensing 
surfaces, 


No experimental procedures were employed other than those outlined 


in the vreceding section, 


Results, 
Table (1), Appendix 2 is a compilation of the above measurements. 
Figure (4) shows typical temperature and yield curves for Still 
FA", plus the water, glass and ambient temperatures, The values for 
Still "B" are indicated where they differ from those of "A", but in 


interest of clarity, the curves are not drawn, 


Interpretation, 

(1) The East and West vertical faces make important contributions 
to the yield of these stills, due to a L/W ratio of one, but in a pro- 
duction=model still the L/W ratio is so large (25-100) that yield of 
the ends is insignificant, 

(2) The South face, since it absorbs 1-2 percent of the incoming 
radiation, is usually several degrees warmer than the North face, and 


consequently could be expected to be a poorer condenser, 


(3) The overall yield coeff, (Yield A+ Yield B) can be taken as 


one when overating under similar conditions, 
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(4) The drop to a coeff. of atout .97 on 22 March and 19 April 
is assumed to be an indication that the inequality due to the afore- 
mentioned tree shadow was not completely neutralized, 

(5) The variation of effectiveness of the condensing surfaces 
and the variation of the yield coeff, on the various days listed is 
also influenced by wind direction, Winds were generally variable in 


the morning, and increasingly westerly in the afternoon, 
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5 Observation to determine the effect on the yield of Still "A" 
of changing the heat transfer conditions of the condenser 


surfaces, 
Discussion, 
tT - T 
The form of the Carnot efficiency equation e = 
Ty 


indicates not only that efficiency increases as4 T increases, but 
also that a drop of X degrees in Tos will increase efficiency more 
than a similar rise in T, since: a0) = Tp 98 

gt, Ty ato 
Likewise with solar energy collectors in general, the sfficiency 
rises as collector temperature decreases, since radiation reception 
is increased and losses are decreased, 

Work by Howe (19) has indicated that the greater portion of any 
energy which may be stored in a solar still at sundown by reason of 
its heat capacity is not regained in the form of yield, so that re- 
duction in the overall teuiperature of a still seems advantageous in 
this respect. 

The enthalpy of saturated vapor decreases slightly as temperature 


decreases ’ 
At 100°F, hfg = 1037.2; hg = 1105.2 


At 150%, hfg = 1008.2; hg = 1126,1 
Howe (19) observed that under conditions of low wind velocities 
and with high incident energy that yield is markedly curtailed, He 
found that in one series of tests the maximum yield with no wind ona 
June day was .05 gallon per day per ft“, while with wind speed above 
nine knots the yield on a similar dey was ,l gallon per day per ee 
Fitzmaurice and Seligman (16) state, "We think that there is 


an appreciable drop in efficiency due to the escape of vapor at the 
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vents, as a corollary of this we assume that evaporation is ahead of 
condensation, and to work at true atmospheric pressure it will be 
necessary to find balance between evaporation and condensation such 
that the pressure buildup is negligible," 

The above facts and observations seem to indicate that advantages 
should incur if the condensing capabilities of solar stills could be 
improved, Such increased capabilities would decrease the overall still 
temperature, and lower the vapor pressure, 

Since such stills are usually oceanside affairs (greatest exception, 
Australia, where brackish well-water is distilled), it seems possible to 
use the sea temperature as a sink rather than ambient air. Although 
this has been proposed in connection with multiple-effect stills using 
focusing collectors, it has not been utilized with stills of the Ias 
Salinas variety to writer's knowledge, 

Although the economic factors involved in the additional pump work 
that would be required have not been explored, it is noted that Bibling 
(26) indicates that high capacity mltiple-effect stills using focus- 
ing collectors should have intermingled flatplate collector stills in 
order to fully utilize the "sun space™®, In an installation of tiuis type, 
the cooling water for the Ias Salinas still could provide feed for the 


nultiple-effect stills. 


Procedure, 

To test the above idea, that is, to increase the condensing 
capabilities, Still "A" was modified as previously noted to allow cool- 
ing the North and/or South condensing surfaces with water, and compari- 


sons were made with the unaltered still, 
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Water at a temperature of 50° - 60°F was allowed to flow over the 
external condensing surfaces, The rate of coolant flow varied on 
various days, but no quantitative measurements were made of the effect 
of this variation, Iikewise, although the temperature of the coolant 
"to" and "from" was measured, no quantitative use of these measurements 


was made, 


Results. 

Although the overall temperature of the still was lowered, and 
although the condensing capabilities were increased, the results were 
in all cases negative, in that the resulting yield coeffcients were 
lower than when "A" was not cooled, Table 2, Appendix 2 indicates the 
compiled results of these tests. Figures (5) - (7) show representative 
temperature and yield relationships for these tests, Only in the first 
hours of a test run, when the temperature of the water and the glass 
surfaces of a normal still are too close together to provide condenging 


potential, did artificial cooling produce any increased yield, 


Interpretation, 

In establishing the reason for failure of external cooling to pro- 
duce improved yield, it might be constructive to establish which factors 
cannot be causitive. 

For this purpose, an analysis is made of the possible losses at 
1200 on 11 April, when both North and South faces of Still "A" were 
cooled, 


Conditions at that time were as follows: 
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still "a" oti 1) ee 


Tt 127% 144.5 oF 

Tn 91 115 

Ts 95 aay” 

Te 98.5 112 

Tw 98.5 112 

Tg, average 94 115 

h 1.738 2055 

? ambient Sine 78.5 

RA + &B 345 
Emissivity, glass and water 296 


Radiation, 
Radiation between collector and glass furface Still "A*® 


L = O7 (Th -T' 
eh 6 Pe ba 28 (i, 


= .939 (.1719 10 ) (se7"- 554") = 38.6 BU, 
HrFt 


Radiation between collector and glass surface of Still "B*" 


similarly: 


g = .939 (.1713 x 1078)" (604.57 - 575+) = 39.4 BIU 
A Erf t© 


Assuming approximate insolation of 200 BLU/Hr ft“ this becomes 


28 or roughly .005 percent increase in efficiency for the cooled 
200 


still over the uncooled, 

While the above value tends to shift to the favor of the uncooled 
still when the sun is further from the meridian, the shift is small, 
and the above figures seem to indicate that radiation is not the cause 


for the discrepancy in yield, 
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Conduction, 
Conduction losses, Still "A" (Considering conduction from tray 


through bottom of still only) 


a = K(ft - Ta) K of fiberglass = .024 
A 58 X = 24" = .208 feet 
= ,024(127 - 78.5) = 5.6 BTU/hr £t2 
208 
Still "B" 
g@ = ,024(144,5 = 78, = 7,62 BIU/Hr rt 
A 208 


or an approximate increase in efficiency of_2 or one percent for 
200 


the cooled still over the uncooled, 


Absorption, 

It was assumed that the absorption of insolation by the thin layer 
of water would not be excessive, although absorption spectrum for 
water rises steeply in the infraw~red region, Overall absorption 
coefficients for the absorption of the sun's energy by water are not 
readily available. Sverdrup (1) gives an overall extinction coeffi- 
cient for seamwater at the surface of 944+, Since by the definition 
of extinction coefficient, this is equal to the absorption coefficient 
at the surface this value was used, Sverdrup also indicates that the 
value for pure sea-water is not markedly different from that of fresh 
water. 

The maximum thickness of the water film did not exceed 2mm, We 


therefore have: 
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~ 944 . on 6944) (002) _ 09981, 


ire = Tod I, = Io 


which indicates that absorption is not an important factor. 


Reflection, 

No quantitative estimate of roflection losses was made, The 
water film, it is believed, did contribute to extra refloction losses 
due to the rippling of the film, particularly on the lower half of the 
glass surface, The water film in this area caused a discernable 
shadow to be cast upon the water in the evaporating tray. However, 
particularly since at most only 40 percent of the glass surface was 
covered with the cooling water film, this extra loss is not considered 
to be the predominate reason for the ,845 yield coefficient, 

Therefore, assuming that the above items are not the dominating 
factors in the reduction of the yield coefficient, the convection pro- 
cess must be examined more fully. 

That convection losses from the water to the condensing surfaces 
in the cooler still must be greater than in the warmer is evident upon 


examination of the mixture ratios of the air-vapor mixture within the 


stills: 
Still "a" Te = Tt +Teg = 127 +94 2 110,5 
2 2 
Pvapor = 1,293 (27) 
Pair = 14,.696-1,.293 = 13.404 


prapor ~ 1.293 x 18 4. » 95975 
Fer 13.403 = 29 





Still "Bf T= +Te = 144.5 +115 = 129.75 
2 2 
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Pvapor = 2.209 (27) 


‘Pair = 14,696 ~ 2.209 = 12,489 


vapor . 2.202 318 = 1097 
p air 12.489 x 29 


Mixture ratio "A" = 05975 = ote 
Mixture ratio "3B" 1097 


The gross assumption mst be that twice as much air must be 
circulated in Still “A® in order to condense an equal amount of water, 
with attendant greater convection losses due to sensible heat transfer 
from the air mass, 

In the extreme cases, of course, mixture ratio of zero would mean 
maximum convection losses, while with a mixture ratio of infinity 
convection losses would be minimum, The cases in between, however, 
are not so patent and will be examined more carefully in the next 


section of this paper, 
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6, Observation to determine the relation of the yleld of the 
condensing surfaces to the dependent variables by traywater 
temperature and external tenPerature of the glass condensing 
surfaces, 

A solar still involves simultaneous evaporation and condensation 
of @ vapor in a mixture of vapor and a noncondensable gas, together 
with convection, radiation and conduction heat transfer, Although each 
field has been subject to mich investigation, the combination, occuring 
in a closed chamber of the geometric shape of a las Salinas Still, to 
the writer's knowledge, has not been so investigated, 

In an effort to establish the relationship between the yield of a 
solar still and the condition existing in the still, it was decided not 
to attempt force-fitting solutions to above individual proglems upon 
the still, but rather to find such a relationship using easily measur=- 
able variables and not involving trial and error methods particularly 
found in existant solutions of the condensation problen, 

The main factors influencing this decision were: 

a) The condensation phase of the problem definitely involves 
both film and dropwise condensation, Any analytical solution would 
have to include a subjective measurement of the relative influence of 
both, 

bd) The convection process within the still is complicated both by 
geometry and because of the dual "potential" of a heated surface and a 
cooled surface within an enclosed space. The convection currents do 
not fully fit any of the standard flat plate free convection assumptions, 
In this connection physical observations of the convection pattern with- 
in golar stills are extremely interesting and easily made, 


In the late afternoon, between 1630 and sunset, if an observer looks 
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into a still from the Hast side with his eye close to the edge of the 
glass side, he will observe a current of tiny water-vapor particles in 
rapid motion moving about the still. These are made visible by the 
reflection of the sun at that hour and are invisible at any other tine, 
So clear is the presentation that these particles can actually be observ-~ 
ed forming, or evaporating, as they move toward a warmer region of the 
still. The convection patterns observed are indicated in Figure (8), 
It is emphasized that this current is rapid, and bears no resemblance 
to the slow rise of vapor from heated pan of water on a cold day. 

The decision referred to above having been made, a standard dimen- 


sional analysis was made, The following variables were considered: 
a,m, kK Ape, €& Co, BT, 9,b, @ DB, A ZB. 
The following nondimensional parameters were obtained: 


Nussult Condensation Nunber 


Nussult Number 


Iie? 
awr_K At 
ia 
K 
—riee Peclet Number 
@ 
PLIA 
PA 
At, Dy 
a’ 
“a ABs 
_e_ 


PHD Schmidt Nomber 
3 
got 

4 n Prandtl Number 


2) 





mys 


The parametersox , € , while important in the design of a solar 
still, are not variable within a certain still, therefore these para- 
meters are discarded, 

Although the condensation number finds its application in the pro- 
blem of condensation of pure vapor, it was thought that the relation- 
ship might be extended to include a gas-vapor mixture. Therefore, an 
attempt was made to correlate this parameter to the Peclet number, All 
properties were evaluated for water at a film temperature which was the 
average of Tt and Tg. The AT used was Tt - Teg. 

It was found that no correlation existed, and this should have been 
expected since the important properties of vapor pressure and mass 
diffusivity are not included, 

A search among other standard heat transfer parameters failed to 
locate one which seems to fully fit the conditions existing in a solar 
still. Conspicuously absent are those containing a pressure tern, 

It was considered that such a parameter must reflect temperature 
and pressure conditions and account for diffusivities. To that end the 
ao combined with the condensation number to obtain 

% 


ahLfPAUe __., This parameter was then separated into the following 
KP oT 


dimensionless parameters: A or” Se oe 
Ce AT x s 


These terms can be considered as follows (7): 


ae A = TIatent Heat 
CeAT Sensible Heat Exchange 


II, ple Dy (lewis #)7> = Mass Diffusivity 
i Thermal Diffusivity 
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ii. gl fP = Gravitational Forces 


P Pressure Forces 


Since mass transfer is perpendicular to gravity flow, the mass 
transfer then should be inversely proportional to the ratio III as it 
is written, therefore, to obtain direct proportionality III was inverted. 
Similarly, mass transfer should be inversely proportional to I as written, 


so IJ was also inverted, 


The parameter then becomes CgAT . Ce . P= 4X, 
A hd yp 

The value of Lewis number for saturated air vapor mixtures has not 
been fully established, but reported values in the region 60% to 1509 
vary from .81 to .91, The variation with temperature over this range is 
slight and confined to the third decimal place, Therefore a constant 
value of .866, reported by Hilpert was used, 

In initial efforts to obtain correlation with this parameter, P and 


p were evaluated for saturated steam at Tf, The pressure parameter 


then becomes RE. No correlation with Pe could be obtained, 
gL 


As finally employed, all properties of I and III are evaluated as 
air-vapor mixture properties at Tf and I] is evaluated as _1. Pressure 
was interpreted as AP = P, - Pe ; 

The* characteristic length was the length of the North or South glass 
condensing surface (47 inches). The parameter Nis most easily computed 
in the form KWEN' AT AP where N' contains all the constants and the temper- 
ature variable fluid properties, 

Pe was similarly evaluated for mixture properties and becomes 4,15 
x ft, 


In Figure (9) are plotted the results of final correlation attempts. 
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While only runs for two different days have been plotted, it is consider= 
ed that due to the varied nature of the conditions imposed, that these 
results show definite correlation, Since correlation was obtained in 

the final hours of this investigation, the writer was unable to examine 
dates other than those plotted. It is fully believed however that data 
on Other dates would correlate, On this plot each of the North and South 
condensing surfaces was separately considered, They will henceforth be 
termed An, As, Bn, Bs. 

These two dates were chosen since on these two dates it is consider= 
ed that the best temperature information was obtained (the ground looping 
of thermocouples on prior dates was previously mentioned), On 29 April 
both stills operated under similar nornal conditions, The sky was com= 
pletely overcast with exception of the period 1530 ~ 1700. 

On 13 June the North face of "A" was cooled, and reflects a corre=- 
spondingly increased yield, while the uncooled South face of "A® reflects 
@ decreased yield, 

In Figure (10) the average Pe and N parameters of the North plus 
south faces are plotted, This plot still does not account for the full 
M of the stills, since the Hast and West faces are neglected. Since Pe 
presents yield in the form of m per foot of length, it was not considered 
advisable to include the triangular ends, 

No plots of Still "B" for 29 April could be included due to lack of 
data for Bs, 

Data and computation sheets for these dates are included in Appendix 
3 Also included are plots of Vapor Pressure vs, Temperature in milli- 
volts and for N' vs, Temperature in millivolts, which were used in calcula- 


ting parameter values, 
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Comments, 

1) The amount of correlation obtained was gratifying to the writer 
Since here an attemvt was being made to correlate data of an averaged 
nature (the mass flow rate) with data of a more instantaneous nature (the 
temperatures), 

Tg is considerably affected by the wind, and in the afternoon when 
winds were usually stronger and gusty, the scatter caused by the sudden 
change in Tg was particularly noticeable, 

2) The temperature information, even on the dates considered, show- 
ed some evidence of being faulty on Still "BE", There were indications 
that the "loop" appearance of the plot of Bn and Bg might be due to 
faulty temperature readings in the afternoon, This possible trouble 
manifested itself in temperature readings, indicating that two adjacent 
terminal boards were three degrees different in temperature, a difference 


difficult to reconcile with observation of the boards, 


te ee ee EE Pe Oe Be Gees: 2k Y WTO Ww . OTT. ce arene 
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Examination of the Parameter N, 
Assuming at this point that there is a correlation between Pe and 


N, an examination of the factors causing the variation of N is indicated, 


Still "B", 13 May, North Side 


Time (1) 0930 PDST (2) 1330 PDST Increase factor 
Col, 2+ Col. l 


Pe 0 7 OF 4 N65 5 ef 
N 77707 2793 15.7 
Tw 94.17 145.23 
Tg 78 110 
c 86 127 5 
QT 16.17 BRS) 2,08 
QP 0330 1.997 6.05 
N! 3330 4010 1.20 


Notes N=wN' pat OP 
Referring to the table above, we see that the most important factor 
in the rise of a solar still to maximum midday production, in terms of the 
variables considered, is the vapor pressure. Vapor pressure is however 
completely expressible in terms of temperature, If P were directly 
proportional to T then the increase in the P term would be the same as 


the AT increase factor, 2.08. The "extra factor", fa08 = 2.91, can be 
2.0 


attributed to the temperature variation, Grouping the temperature factors 
we then have 2.91 x 1.2 = 3.5. Similarly the 4? factors = (2.08)* = 
4.325, meaning that variation inOT is roughly a 20 percent greater factor 
on this particular day, However, on a hot day with low wind, the tempera- 
ture factor would tend to increase, and the AT factor decrease, 


Perhaps & more meaningful examination would be to examine the change 
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in O T necessary to offset a change in fT, This requires a trial and 
error approach since 4 P is affected both by 47 and Tr, and it is 
found that at 1330 if the film temperature is dropped ten decrees F,. 
then 4 T must increase 4,8° to a total of 40°F, to maintain the same 
yield. In terms of the water temperatures and glass temperatures, if 


Tt drops ten degrees, Tg must drop 16,3° to maintain the same yield, 


Interpretation of N — Pe correlation, 

While the N - Pe correlation might be more complete if an interior 
glass temperature were used, this was avoided in an effort to finda 
correlation using the most easily measured important dependent variables, 
Since an exterior glass temperature was used it can be expected that if 
a thinner glass, or one with higher thermal conductivity were used, that 
the N — Pe curve would have a slightly greater slope. 

It is believed that the N - Pe correlation is independent of 
efficiency, i. e., for a certain combination of water temperature and 

QT the mass flow rate will be determined, 

It is hoped that the value of the above correlation lies in four 
regions which will be separately discussed, 

1) In the method of variation of the parameter itself: 

The fact that OT must increase sharply to offset a decrease in 
the water temperatures indicates the probable reason for the failure 
of the efforts to increase efficiency by increasing the capability of 
the condensing surfaces, While no method was found to measure the con- 
vection losses, due to inability to measure the insolation, it is con- 
gidered a reasonable assumption that they vary as the 5/4 power of A T, 


in some close accordance with the usual formas for the free convection 
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of air under a cooled flat plate, where h = -27(a 2) 02) With this 
assumption, we see that any attempt to cool a solar still is likely to 
cause increased convection losses, since the overall effect of that 
cooling is to decrease the water temperature and increase the AT, We 
further see that such cooling would increase the efficiency of a still 
only in the operating region where the increase in convection losses is 
offset by a larger decrease in radiation and conduction losses, This 
probably occurs under some combinations of high ambient temperature, 
high insolation, and low wind velocities, It should seldom occur on 
the Monterey Peninsula, where the ambient temperatures rarely reach 70°F, 
Assuming the proportionality of our convection losses to the 5/4 power 
of ST, we see in the preceding example where the Tt was lowered ten 
degrees requiring the T, to decrease 16,3°, that the ratio of the con- 


vection losses of the final and initial conditions becomes 


(ita = 242) = 1.23, or a 23 percent increase in convection losses, 

While it would be hazardous to say that the N - Pe correlation under 
discussion can be applied to a still using forced convection, probably 
gome very similar correlation can be obtained, Daniels (4) outlines an 
experiment with a forced=-convection type of still and concludes that the 
low efficiency (25 percent) was due to too rapid air movement which re-~ 
sulted in the air not being saturated. This writer is in agreement that 
the air circulation may have been too rapid, but suggests that perhaps 
the result was not unsaturated air, but rather a lowering of the system 
temperature, together with the attendant increased f T and convection 
losses, 

Similarly noted in this same gource (4) is a plastic lifeboat still 


developed by Telkes, which floats on the ocean surface, It is suggested 
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that this still might produce greater yield if it were lifted out of 
the water, 

Fitzmaurice (16) indicates an increase of yield of 44 percent when 
& still with a diamond=shape cross-section, with the tray suspended in 
the center, was used instead of a conventional Ias Salinas type. Al- 
though some of the increased yield no doubt is due to the lack of con= 
duction losses, he attributes this increase to the enlarged condensing 
capabilities, and to check tnis assumption, he covered the lower surfaces 
of the diamond with insulating board, whereupon he noticed a "small but 
definite reduction in yield," 

Examination of the included plot of brine temperature, ambient 
temperature and yield curves for the las Salinas type gives strong 
indication that these stills were in fact operating under conditions 
where a reduction of the water temperature by increased condensing 
capability should increase efficiency, Ambient temperatures for these 
measurements were often above 100° F, It can be questioned, however, 
whether the diamond=shape stills in seasons other than the summer 
season reported on, will show increased efficiency beyond that due to 
reduced conduction losses alone, 

As a corollary it would seem that an effort to increase the operat= 
ing temperatures of a solar still should pay dividends, up to the point 
where the rising conduction and radiation losses overcome the effect of 
decreased convection losses, 

If the above surmizes are correct, perhaps the most important con- 
clusion to be drawn is that various operating locales require different 
still designs, and that perhaps there ig no one best design for all 


locales and seasons, 
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2) The last paragraph serves to introduce the second major region 
where the N parameter may prove useful: 

It is intuitively believed that the performance of a solar still is 
reflected in the relationship of the temperature of the water to theA T 
vetween water and condensing surface, particularly under conditions des- 
eribed by Howe (19) of low wind and high insolation, where the performance 
radically changes, Therefore, it is suggested that if other investigators 
can confirm, improve, and extend the correlation described here, a whole 
new method of solar still design and operation might be evolved about 
the parameter N which will allow designing a still for specific locales 
and ambient conditions, 

Since N is a function of T and AT only, using the coordinates of 

OT and T, the family of constant N lines can be drawn, 

Since 2 line of constant Nis similarly a line of constant m, then 
the N - Pe correlation permits relating in terms of T and AT, the actual 
instantaneous performance of @ solar still to the theoretical maximum 
performance, That is, the actual operating point may be studied in rela- 
tion to the theoretical maximum N curve, 

It is considered possible, particularly where means are available 
to measure insolation, that further study might suggest a method of 
laying out upon these constant N curves the family of "mximun efficiency 
operating curves" for a solar still, As the film temperature and DT 
varied throughout the day, departure from these curves might indicate 
design or operating changes to bring about closer adherence to one of 
the curves of tnis family. The gradients to the family of N curves 
might be found to have some significance in this respect. 

It is regretted that time is insufficient for this writer to pursue 
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this idea, but it is hoped its interest to other investigators might 
be such that it will in the futuré be proved or disproved, 

3) The third region where the parameter N may prove of value is 
as a substitute for m in any correlation where m might be employed: 

For example if efficiency is under consideration a study of N should 
prove as effective as a study of ad, and have the great sina ttcestttt 
all quantities are instantly measurable, do not require consideration of 
prior or future measurements, as any determination of m generally does, 
and time lag within the system does not effect the correlation as it 
generally does with correlation of mass flow rate, 

4) The fourth region where N may have value is in other related 
fields of heat transfer, 

It is considered that not only can finer Agasuseneath, together 
with an inclusion of a "time lag" factor increase the correlation of NWN 
with Pe for a solar still, but it is expected that this parameter may 
find applications in other fields of heat transfer theory where masa 
transfer of condensable vapor from a noncondensable gas is considered, 

While the writer is not sufficiently familiar with heat~-transfer 
literature to say with assurance whether this parameter or one similar 
to it has appeared previously, he notes the comments of Klinkberg & 
Mooy (7) on the pressure parameter of N, (_P_) to the effect that in 
@ll analogy theories published so far, Pe mentioned ratio is dis=- 
regarded, While there may be excellent reasons unknown to this writer 
for tais, it is hoped that this parameter has now found at least one 


useful application, 
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APPENDIX I 


EXTERNAL ENERGY BALANCE FOR STILL "A*® 


at 1200 on 19 April, 1958 





DEVELOPMENT OF ONE SLANT FACE AND 1/2 OF EACH VERTICAL END FACE, 


Tt = 139° F, eens 2.5 knots, 090° 
een 112 3", | = 2.268 lbs/hr. 
T(dry bulb) = 64°F, A = 1014 BIU/}be 

T(wet bulb) = 60° F, Ag = 28.24 £t*, 


Ag less E and W faces = 21,52 ft” ° 


Energy removed from still by wind: 


Wa = W (12) 
(Kair(Tg-Ta) 


Nu = .592 ¥ Re (2) 


Re = V LP 


Aa 


L(effective) = C + 2(2/3D) 


= 47/12 + 2(2/3) sa 5.642 ft. 


Effective width, W= Ac = 28,24 f¢° = 5.01 ft. 
LT 5,642 £t 
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Nu = 2592 \/(205 : 1,689) ° 5.642 = 231 
156 10 


q/W = 231 Ka(Tg - Ta) 
= 231(,016)(112 - 64) = 177.3 


q = 888 BTU/Hr, 


Energy removed from still by natural convection? 


For a vertical plate: h, = .29 (ap) *? (2) 
L 2 


For a heated horizontal plate: h, = .27Q2)° 
L 


Uses h, = 26 (0R)"*? for 45° plate, 


LN&S = Slant height of face = 33 feet = 2.75 feet. 


LE&W = 2/3 altitude = 2/3 22/12 = 1.22 feet. 


L average (area weighted) = 2,38 feet 


ho = .280T)°*9 = 2255 b T°~9 
2.38 

q/Ag =he P= ,2255(Te - Ta)" 

= 2255 « (4B)t29 


= 28.64 


q = 810 BIU/ E 


Energy removed from still by radiation to sky? 


For N and S faces only; assume zero for BE & W faces, 


by 
qg/A= Tabs(.22 + .148 ° “iui eb) (5) 


Tw/Ta = 60/64 .°. Rel. humidity = 80 percent 


oe Pw .24 PSIA = 12.44 MM Mercury 
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Tg mw 113° FF. = 45° C = 318° K 


qQ/A = ,.826 ° 1074 (318°) (.22 + 148 © 1972068 © 12M) a5 


= .1735 = 36.4 BIU/Erft’. 


A = 21,52 £t* q = 827 BIU/Hr. 


Heat removed from still by conduction: 


Bottom of still only, considered, 
Bottom insulation, 2k" fiberglass blanket, assume K = ,024 


q/A = K(Tt - Ta) = 02% * (75) = 8.65 BIU/ft“Hr. 
x Lad 
12 


A Tray = 15.34 ft 


Summation: 


q wind = 888 BIU/Hr 
q convection = 810 BIU/Hr 
q radiation # 827 BIU/Hr 
q conduction = 133 BIU/Er 


2658 BrU/Hr 


Por sun at midday, reflection losses = 8 percent of incoming 
insolation, 


Q total = 2658 = 2890 BrU/Er 


092 
Collector area = 15.34 ft* 
Q, per unit area of collector = 188.2 BrU/ft“Hr. 
A = 2,268 © 1014 =» 2300 BIU/Er 
= 150 BRU/Hr per unit collector area, 


Midday efficiency = mg=150 = 9.0 percent 
Q 188.2 a ee oe 
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APPENDIX 2 
Compilation of Results for 
Sections 4 and § 
Key to Tables 1 and 2: 
Column 1 Still identity. 


2 Amount.of water condensed per unit area by North face, 
2 
los /ft~. 


Area North face each still = 10.78 ft’, 


3 Amount of water condensed per unit area by the South 
face, lbs/ft*. 
Area of South face of each still # 10.78 ft’, 


4% Amount of water condenged per unit area by the Hast 
plus West faces lbs/ft*, 


Area of East plus West faces = 6,72 ft’, 
5 Total yield per unit area of collector surface, 1bs/ft*. 
Area of collector = 15,34 ft”, 


6 Coefficient of performance of Still "A": 
total yield "A® 4 total yield "B", 


7 Comments. Times indicate cooling period, 
OC indicates overcast, 


val 





Date 


31 January 


4 February 


5 February 


6 February 


9 February 


10 February 


ll February 


9 March 


ll March 


18 March 


22 March 


19 April 


td & 


TABLE I 


Compilation of Results for 


section 4 


5 


oO5 
035 


» 264 
0264 


2267 
0269 
0147 
0143 


563 
oD D0 


0392 
e351 


0354 
0354 


0663 
2670 


0652 
0647 


soe 
0635 
0962 
580 


2910 
0934 


Average Coeff of Performance: 


42 


1,03 


101 


1.003 


1,00 


991 


1,006 


0995 


0970 


972 


1.004 


7 


Thin overcast 


90° overcast 


Clear until 
1000, then 
overcast 


Thick over= 
cast 


Broken clouds 
Overcast until 
1200, then 


clear 


Broken over@- 
cast 


Clear day 

Very thin over- 
cas t 

Thin overcast, 
thickening at 
1415 


Thin haze 


Clear, ho 








Table 2 


Compilation of Results for 


section 5 
Date 1 2 3 Ly 5 6 7 
North Face of Still A Cooled 
26 February A 9482 241 e2l9 9003 0951 1330-1600, over= 
B o422 2314 ~206 2034 cast after 1500 
7 March A 428 0220 0220 SO 7 0970 0930-1700, over= 
B 2306 2265 olle2 o 584 cast after 1430 
19 March A 2493 e208 e210 0002 e945 0915-1700, medium 
B 0 389 EV SASS. 0035 overcast in P.M, 
be April A ~705 eO05 0254 819 940 0905-1700, clear 
B sO 7 0453 0433 872 Ta max, = 81° 
13 May A 6765 437 0355 1,000 0929 0915-1845, max. 
B Aon 588 0525 1.077 cooling flow - 
thin haze, clear 
Average Coefficient 0947 sahara 
Soutn Face of Still A Cooled 
5 March A ct 2260 Paha 0434 0910 0900-1530, over- 
B 0297 e2tt2 gees 477 cast after 1500 
10 april A RS ies, “51e0) 0296 ata é 6909 0930-1700, over= 
B ABs 3 443 2427 860 cast after 1700 
Average Coefficient 907 
North and South Faces of Still A Cooled 
6 March A 5336 2.326 «vl80 0 54 869 0900-1530, over 
B 2595)! (Sele 2.95 2025 cast after 1500 
mieapril A 5120 61 >» 187 e775 .880 0930-1705, clear 
B So: 452 441 sO7o Ta max, = 75° 
Average Coefficient a07 5 
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APPENDIX 3 


DATA SHEETS FOR 


29 April and 13 May 1958 
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